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In fully developed three dimensional fluid turbulence the fluctuating energy is 
supplied at large scales, cascades through intermediate scales, and dissipates at 
small scales. It is the hallmark of turbulence that for intermediate scales, in the 
so called inertial range, the average energy flux is constant and independent of 
viscosity pj-lsl]. One very important question is how this range is altered, when an 
additional agent that can also transport energy is added to the fluid. Long-chain 
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ymers dissolved at very small concentrations in the fluid are such an agent 



7| we introduce a theory 



51]. Based on prior work by de Gennes and Tabor 
that balances the energy flux through the turbulent cascade with that of the 
energy flux into the elastic degrees of freedom of the dilute long-chain polymer 
solution. We propose a reflned elastic length scale, r^, which describes the effect 
of polymer elasticity on the turbulence energy cascade. Our experimental results 
agree excellently with this new energy flux balance theory. 

To date, the only relevant theory on the interaction between polymers and turbulence 
cascade is the "energy balance theory" by de Gennes and Tabor {g], Q] as summarized recently 
by Sreenivasan and White js]. It states that the turbulence energy cascade is essentially 
unaltered down to a scale at which the energy stored in the elastic degrees of freedom of the 
polymer is equal to the kinetic energy of the flow. So far experiments on fully developed 
three dimensional turbulence do not convincingly support this theory 9|. Here we argue 
that this may not be surprising as it is the turbulent flux of energy and not the energy itself 
that determines the inertial range properties of turbulence. 

In this paper, we provide what may be called an "energy flux balance theory". In our 
theory, the turbulent energy flux through the cascade (or the turbulent energy transfer 
rate from scale to scale) is gradually reduced by the energy transfer through stretching and 
recoiling of the polymer chains, with the elastic energy flux becoming dominant at small 
scales. For the sake of making progress, we assume that the balance of elastic and turbulent 
energy occurs in average at one length scale. This of course is a crude assumption that does 
not hold in detail as the turbulent energy flux is known be intermittent js]. Nevertheless, 
we may expect that it captures the main features, just as Kolomogorov's 1941 (K41) theory 
does for pure fluids. Here we show, that the energy flux balance theory is supported 
quantitatively by our experimental data measured over a wide range of parameters in fully 
developed turbulence. Our approach may also be applied to and provide new perspectives 



into turbulence in other complex flows involving two or more nonlinear mechanisms, like for 
example, turbulence in conducting fluids, plasmas or quantum fluids. 

For stationary turbulence in Newtonian fluids, the kinetic energy of the turbulent motion 
per unit mass is supplied at a rate of Sj at the forcing scales: length scale L and time 
scale Tl, which are related by ej ~ /T^ ~ u^/L, where u is the root mean square of 
the fluctuating velocity. As the scale decreases, the viscous effects become more and more 
important. It is well known, that at the Kolmogorov scale the viscous forces dominate 
and dissipate kinetic energy to heat at a rate of Sn^ which gives naturally the length scale 
r] = {u^/sdY^^ ^he time scale = {u/soY^'^, where u is the kinematic viscosity of the 
fluid. The Reynolds number Re = uL/u, a parameter to characterize the intensity of the 
turbulence, measures the separation of scales: L/r] ~ Re^^^ and T^/Tri ~ Re^^^. For Re ^ 1, 
the forcing scales (L and Tl) are widely separated from the dissipative scales {rj and r^). 
The K41 theory [2| then states that for intermediate scales 77 <^ r <^ L or ^ r ^ T^, 
the energy is transferred down to smaller scales without loss. This immediately leads to the 
conclusion that ej = eT{r) = en = e is a. constant. Here we used the notation exir) to 
emphasize that the energy transfer rate depends on scale r. This local energy transfer rate 
can be estimated as erir) = Ur'^/Tr, where Ur and are the characteristic velocity and time 
at scale r {ur is often related to the velocity differences at scale r: Ur ~ \u{x + r) — u{x)\). 
The results from K41 theory is Ur ~ {erY^^ and Tr ~ which is consistent with 

^rif') = £ for 7] < r < L. 

This elegant picture of the energy cascade is changed when flexible long-chain polymers 
are dissolved in the fluid. For simplicity, one may regard a single polymer chain as an entropic 



spring that is constantly stretching and coiling back in the flow [6|, 17|, llOU . The turbulence 
fluctuations at different scales contribute unequally to the stretching of the polymer chain. In 
particular, Lumley concluded that only those fluctuations with time scale < can stretch 



the polymer chain [llj, where Tp is the entropic viscous relaxation time of the polymer chain. 
This "time criterion", which is essentially the same as requiring that the scale-dependent 
Weissenberg number Wi^ = Tp/r^ > 1, thus deflnes the Lumley scale r* = {eXp'^Y^'^ (see e.g. 
[sl). The physical meaning of r* is that below this scale the local fluid deformation would 
be strong enough to stretch polymers. 

Please note, that the Lumley scale can only tell us on whether a single polymer chain 
can be stretched by turbulence or not and thus cannot address how polymers dissolved at 



a certain concentration will affect the ffow. This was addressed by the "energy balance 
theory" proposed by Tabor & de Gennes. By analogy to polymers in a linearly stretching 
field (see Methods), Tabor & de Gennes {g] suggested that the polymer elastic energy per 
volume is Ef,{r) ~ CpkT{r* /rY"^"^ for r <^ r*, where Cp is the number of polymer chains per 
unit volume, k is the Boltzmann constant, T is the temperature of the fluid, and n is an 
unknown exponent that is related to the average stretching dimensions of the local flow field. 
De Gennes 7| further argued that the turbulent energy cascade will be truncated below a 
scale r** at which the polymer elastic energy balances the kinetic energy of the turbulent 
fluctuations (see also Figure [T|^a)): 



pul„ = p{eTT**f" = Cpkn^f-'\ (1) 



This gives 



r** = {kTp-^CpeT~~^Tp—)J+^. (2) 



The polymer relaxation time ^ (N^^^aYp/kT, where is the number of monomers per 
chain, a is the length of a monomer, and p is the dynamic viscosity of the fluid. 

When writing down Eq. ([1]), de Gennes conjectured that the turbulence energy cascade 
is unaffected at scales r > r**, which, however, is not consistent with the theory itself as 
it assumes that polymers already gain elastic energy from the stretching by eddies of size 
r** < r < r* and hence must have diverted part of the turbulence energy flux at scales 
r > r**. Assuming the time scale for polymers to transfer elastic energy down scales is Tp, 
then the elastic energy flux is 

JcTr r* 

.e(r) ~ E^{r)/irpp) ~ ^(^)^"/^ (3) 

TpP T 

The elastic energy flux given above increases as r decreases and can dominate the turbulence 
energy flux Et, as suggested in Fig. [T](b). A new scale can be defined when the two fluxes 
balance: 

eT = e,{r,)=A^C-r^', (4) 

TpP Te 

which yields 

/ timl \ -2- -2- 1 2_ 3__2_ . , 

= [AkT / p)^^Cp^'-^Srp'^ ^"Tp2 

where a proportional factor A is introduced to convert the scaling relation to an equation. 
For < r < r*, the turbulent energy flux is only slightly affected and inertial range scaling 



of the cascade would still capture the behavior to leading order. However, for r < r^, the 
turbulent energy transfer will be strongly reduced. 

As a test of this conjecture we investigated experimentally the validity of the elastic energy 
flux as given by Eq. |H The experiment was done in the von Karman swirling (VKS) flow 
system (see Methods for details). We probed the energy transfer in the flow by measuring the 
second order transverse velocity structure function D^fyflr) = (|n_L(x + r) — n_L(x)p), where 
u± is the velocity component in the direction perpendicular to the separation vector r. For 
Newtonian fluids, in the inertial range, D]\ij\f{r) = |C2(£t'")^^'^; where C2 is a constant (here 
we use C2 = 2.11). This expression can be rewritten as D]\fN{r) = ^[jDj^^{r) /C2]'^^'^ = St- 
Figure [2] shows DN^^r) for the flows at fixed Rx = 360 with different polymer concentrations. 
{Rx is a Reynolds number definition widely used in the turbulence community, which is 
related to the more familiar Reynolds number as Rx ~ Re^^'^.) For the pure water case, 
DNNij") displays a plateau of value St in the inertial range. When polymers were added, 
DNN{f) reaches a plateau at larger r and the value of the plateau is lower. Both the 



suppression of small scales ITj and the decreasing of Et |9| have been observed before. 

In |9|, it was also noticed that et measured from inertial range exceeds €d measured from the 
dissipation range by a large amount, which is consistent with the conjecture that at small 



iuctuation energy by elasticity and 
9| a critical polymer concentration 



scales it is the polymers that transfer part of the fluid's 
hence reduce the turbulence energy flux. Note that in 
of 0c = 7 ppm was observed below which et was not affected by the presence of polymers 
in the flow with similar Rx in VKSl. In VKS2, Et was found to decrease slightly even 
at = 1 ppm. This difference is most likely due to the change in the large-scale flow, 
as a consequence of the modifications in propeller size and the vane structure. As we will 
show next, when data measured from the two apparatuses are processed the same way, the 
results overlap with each other, which strongly supports the conjecture that the mechanism 
of elastic energy transfer by polymers is independent of how energy in injected into the flow 
at large scales. 



' DNN{r) shown in Fig. |2] gives the energy 
9[, for VKS flows polymers reduce large 



In principle, the difference in the two curves o 
flux by polymers Ee{r). However, as shown in 
scale velocity fluctuations also, which contributes to the decrease in et- To account for this 
change we assume that the plateau values of DN]\r{r) for polymer solutions correspond to 
the total energy transfer rate Et- The difference between Et and DN]\f{r) at smaller r is thus 



ee{r). As we shall see later, this assumption is well supported by the observed collapse of 
the data. Moreover, we notice that Eq. ([3]) can be rearranged as 



kTc 



V J 



A{r*/r) 



(6) 



From our data we fit the constant A and the exponent n using the measured [^^^^^^^ in 
the form of Eq. In the inertial range our data support the power-law relation for scales 
r ^ rj, i.e., [^^^^Y^^ ~ (r/r*)". From the 14 data sets, we obtain n = 1.0 ±0.2. As shown 
in Fig. [31 all the 14 data sets are consistent with (r/r*)~^ scaling for relatively large r. In 
the inset of Fig. [3l we show the values of A fitted with n = 1 from each data set as a function 
of Wi. Except for the two small Wi cases, A is within 101 ± 17 throughout the parameter 
ranges we explored. The larger values of A found for the two low Wi cases {Wi = 1.0 and 
1.9) is very likely due to that the Wi are below a critical Wic recently found in numerical 



simulation where Wic is reported to be 3 ~ 4 

Previously some of us jo] had found that data from varying Reynolds number but constant 
polymer concentration {Rx = 200, 240, 290, and 350, = 5 ppm) could be collapsed with 



the normalization scale ~ Wi~^'^^ [sl, which is equivalent to ~ st"^''^^ ot n = 0.51. The 
range of Wi in this previous work was relatively small {Wi = 1.0, 1.9,3.4, and 5.8)and the 
change in the value of A (see inset of Fig. |3]) masked the analysis. 

According to the conjecture by Tabor & de Gennes, the exponent n = 1.0 implies that 
polymers are most effectively stretched in locally biaxial extensional regions of the flow. 
This is consistent with the properties of velocity gradients in turbulent flows and previous 
findings of polymer behavior in turbulence. In homogeneous and isotropic three-dimensional 
turbulence, locally biaxial stretching is more probable than uniaxial stretching 19[. More- 
over, recent simulations showed that polymer extensions are larger in the biaxial stretching 
regions of a channel flow 20[ or a homogeneous shear flow 2]J]. Furthermore, we note that 



a previous compilation of drag-reduction data from different turbulent pipe flows reported a 
value fti 2/3 |8[ that is quite close to what we found here from velocity structure functions 
in bulk turbulence. 

With n = 1, we can determine the scale from Eq. ([5]): 



A{kT/pr%'-'eT'W-'- 



(7) 



In Fig. m we plot [^^^^^] ^ /[A{r* /ry-^] against r/r^. All the 14 data sets collapse and show a 



plateau when r ^r^:, which is the expected scahng range for elastic energy flux. For r < 
polymer elasticity dominates the energy transfer process, and the traditional K41 energy 
cascade does not hold anymore. Thus at those scales the elastic energy transfer rate could 
not be obtained by simply subtracting Di^iyf{r) for > from the = case. 

In addition to collapsing the elastic energy flux, should also collapse the structure 
function data since is the scale at which polymer elasticity dominates the turbulent 
energy cascade in the inertial range. Figure |5] (a) shows DNN{r), which is defined as -Dattv 
normalized by the corresponding et- Clearly the lower ends of the inertial range are pushed 
to larger scales when polymers are present, because at the smaller scales the polymers are 
stretched and truncate the K41 type energy cascade. When r is normalized by r^, DNj\f{r) 
from different data sets collapse, as shown in Fig. 5(b). We also performed experiments 
with a different polymer - Polyethylene Oxide (PEO, MW = 8 x 10^). The measured -D^rAr 
at Rx = 342 for two different concentrations [dppm and lOppm) are found to collapse with 
the PAM data but with a different value A ^ 25, which is consistent with the theory as the 
numerical factor A depends on the combination of polymer and solvent. 

In a certain sense, is similar to the Kolmogorov scale rj as both correspond to a scale 
at which the inertial range turbulent energy cascade is truncated by a mechanism whose 
effect is negligible in the inertial range but increases at small scales. On the other hand, r) 
for Newtonian flows and for polymer solutions are significantly different. For Newtonian 
flows, the velocity field below rj is smooth and the small scale turbulence is expected to 
be universal. For polymer solutions, at scales r < there is still room for interesting 
dynamics. It has been shown that in a smooth velocity field the elastic instability can drive 
polymer solutions to elastic turbulence Our theory suggests that in turbulence this 
might occur when the polymer stress dominates the fluid stress, i.e., at scales below the 
de Gennes scale r**. Numerical simulations of decaying turbulence with polymers indeed 
showed an enhancement of energy spectra at small scales (smaller than 77) 22|. It will be 
very interesting to study the small scales experimentally with suitable diagnostic techniques. 



Methods 



Elastic energy of polymer in turbulent flow The behavior of a long-chain polymer is 
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similar to that of a spring, except that the restoring force for the polymer chain is the 
entropic force. The elastic energy in one chain is ~ A;T[A(r)]^/^, where kT is the "spring 
constant" of the polymer and A(r) is the extension of the polymer due to eddies of size 
r. The polymer elastic energy is proportional to A^^^ instead of A^ for a harmonic spring 



because of additional repulsions between monomers 23|. The elastic energy per unit volume 



is then ~ Cp/cTA^^^. To formulate the polymer extension in a turbulent flow, de Gennes & 



Tabor 



n rf 



7| considered the polymer extension in a laminar converging flow, where the flow 
velocity and shear rate are l/(47rr'^~^) and l/{2'Kr^) respectively, where D is the spatial 
dimension of the converging flow: D = "i means flow into a capillary tube and D = 2 
means flow into a long slit. The polymer starts to be stretched when the shear rate exceeds 
its relaxation 1/rp. Thus l/{2TTr^) = l/xp gives a scale r*: the polymer will be stretched 
when it is in the range r < r*. The extension of the polymer is A = (r*/r)-^~^ since the 



polymer extends passively (affinely) with the deformation of the local volume element [24 1. 
de Gennes & Tabor then transferred this idea to the case of turbulent flows and assumed 
that the polymer extension can be written as A = (r*/r)", where r is now the eddy size 
and n is an unknown parameter given by the average stretching dimensions of the local flow 
field. 

Experimental setup and techniques The turbulent flows were generated by two counter- 
rotating baffled disks in a cylindrical tank, known as the von Karman swirling flow. Two 
VKS apparatuses were used in the experiments. One has been described in details be- 



fore [25|]. The diameter of the disk was 25 cm. The diameter and the height of the tank 
were 49 cm and 63 cm, respectively. The other VKS had nearly exactly the same but with 
slightly smaller disks (20 cm of radius). The vane inserts used to break the large scale flow 
in the tank were also slightly different. Using Lagrangian particle tracking technique QjH], 
we tracked simultaneously hundreds of tracer particles in a measurement volume of approx- 
imately (3 cm)^ at the center of the flow, where the turbulence is close to homogenous and 
isotropic and the fluctuating velocities are much larger than the mean. The velocities of 
the particles were then obtained by differentiating the particle trajectories. The polymer 
used was polyacrylamide (PAM, molecular weight 18 x 106, Polysciences Inc.). In the ex- 
periments, Rx = (ISu^/ez/)^/^ was varied between 200 and 360, Wi = Tp/r^ from 1.0 to 
13.0, and polymer concentration from to 10 ppm (parts per million by weight). For the 
experiments with polymer solutions, we used the same R\ as that measured from the pure 
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water case (0 = 0) at the same rotating frequency of the propellers. 
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FIG. 1. (a) A sketch showing the balance between turbulent fluctuation energy and the polvmer 

n a 

elastic energy, which gives the scale r** [6|, lZ[; (b) Similar sketch for the balance between the 
turbulent energy flux and the elastic energy flux, which determines the scale r^. 
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FIG. 2. Compensated second order transverse velocity structure functions Dattv at Rx = 360 with 
different polymer concentrations. 
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FIG. 3. The measured elastic energy transfer by polymers Se{r) in the form of [-^7^^]^ vs. 
r/{Ar*) for 14 different data sets. Inset: the numerical proportional factor A obtained from fitting 
each of the 14 data sets to [^^^]^ = A{r/r*)~^-^. 
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FIG. 4. Collapsing of the measured elastic energy transfer by polymers when r is normalized by 
Tg. In this compensated plot, ([^^5^^]^ is divided by A{r* /r)^-^) and r is rescaled by the elastic 
scale r^. The legend is the same as in Fig. [3l 
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